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We present a microscopic model for coherent pion production off nu- 
clei induced by neutrinos. This model is built upon a model for single 
nucleon processes that goes beyond the usual A dominance by including 
non resonant background contributions. We include nuclear medium ef- 
fects: medium corrections to A properties and outgoing pion absortion via 
an optical potential. This results in major modifications to cross sections 
for low energy experiments when compared with phenomenological models 
like Rein-Sehgal's. 
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A proper understanding of neutrino-induced pion production off nuclei 
is very important in the analysis of neutrino oscillation experiments. For 
instance, 7r° production by neutral currents (NC) is the most important v^- 
induced background to — > v e oscillation experiments, [Tj. Similarly, ir + 
production by charged currents (CC) is an important source of background 
in Up —s- u x disappearance searches [2]. We will follow [3] to describe the 
coherent CC pion production reaction induced by neutrinos 

vi(k) + A z \g S {PA) -» r (A/) + A z \ gs (p' A ) + ii + {K) , (1) 

where the nucleus is left in its ground state, in contrast to incoherent reac- 
tions where the nucleus is broken or left on an excited state. 

We build upon a microscopic model for the single nucleon process (vN — ► 
l~Nir + ). We sum coherently the contribution of all nucleons on the initial 
nuclei, which is modeled after a Fermi gas in Local Density Approximation. 
Coherent ir production is most sensitive to the Fourier transform of the 
nuclear density for momentum q — k^, which gets its maximum value when 
q and k n are parallel. For this particular kinematics the vector contribution 
to the single nucleon (W + A — ► Ntt) currents, which is purely transverse 
k w x q, vanishes unlike the axial contribution. This dominance of the axial 
contributions is exploited through the PCAC hypothesis by the Rein-Sehgal 
(RS) model [U E], to relate the neutrino coherent pion production cross 
section with the pion-nucleus elastic differential one. 

For the elementary process we use the model derived in 
Ref. [6], see Fig. [Q In addition to the 
A(1232) pole (AP) (first row) mech- 
anism the model includes background 
terms required by chiral symmetry: nu- 
cleon (second row) pole terms (AP, 
CAP) contact (CP) and pion pole 
(PP) contribution (third row) and pion- 
in-flight (PP) term. Background terms 
turn out to be very important and be- 
cause of them, the flux-averaged u^p — > 
fi~pir + ANL cross section [TOJdl] is de- 
scribed with an axial form factor where 
the dominant C§ nucleon-to-A axial 
form factor was fit to data resulting in 
C^(0) = 0.867 and M AA = 0.985 GeV. 
This value for C^(0) is significantly 
smaller than the value of about 1.2 de- 
duced from the Golberger-Treiman re- 
lation (GTR) used in PCAC-based approaches like RS. 
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Fig. 1. Model for the W+iV -> N'n re- 
action. The circle in the diagrams stands 
for the weak vertex. 
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When applied to a coherent process in finite nuclei we find that the NP 
and CNP nucleon pole term contributions partially cancel each other, that 
the PF term does not contribute to the coherent cross section 
and the CT and PP terms van- 
ish for isospin symmetric nuclei. 
As seen in Fig. [2] the effect of 
the background terms, both in the 
plane wave impulse approximation 
(PWIA) and in the full model cal- 
culation, is very small. Thus, we 
predict cross sections around a fac- 
tor of (1.2/0.9) 2 ~ 2 smaller than 
approaches assuming GTR. In the 
following we will always use the full 
model of Ref. [6] with C^(0) = 
0.867 and M AA = 0.985 GeV. 

Nuclear medium corrections to 
the dominant A diagram are considered by including the self-energy of the 
A in the medium, Ref. f7|. Another major nuclear medium effect is pion dis- 
tortion, which is taken this into account by replacing the plane wave with a 
pion wave function incoming solution of a Klein- Gordon equation with a mi- 
croscopic optical potential, Ref. [8]. In left panel of Fig. [3] we show the pion 
momentum distribution (LAB) for CC coherent pion production, in the peak 
energy region of the T2K experiment. Including A in-medium self-energy 
(long-dashed line) reduces the PWIA results (short-dashed line). Further 
inclusion of pion distortion (full model, solid line) reduces the cross section, 
and the peak is shifted towards lower energies. The total cross section re- 
duction is around 60%. Medium and pion distortion effects in coherent pion 
production were already evaluated in Refs. [9j. However, the authors of these 
references neglected the nucleon momenta in the Dirac spinors. The effect 
of this approximation (nucleons at rest, dotted line) results in a ~ 15% de- 
crease of the total cross section. In the right panel of Fig. [3] we show the pion 
angular LAB distribution with respect to the incoming neutrino direction. 
The reaction is very forward peaked, as expected due to the nucleus form 
factor. The angular distribution profile keeps its forward peaked behavior 
after introduction of nuclear medium effects. 

We examine in Fig. [4] the NC differential cross section with respect to the 
variable E^fl — cosO^), proposed by MiniBooNE. Our prediction is appre- 
ciably narrower than that displayed in Fig. 3b of Ref . Q]. The MiniBooNE 
analysis relies on the RS model, so we try to understand the differences 
between this and our model. RS's expression for the coherent tt° produc- 
tion cross section was deduced in the parallel configuration, for which the 



All diagrams; Full model 
Only A; Full model 
All diagrams; PWIA 
Only A; PWIA 




Fig. 2. Pion momentum differential LAB 
cross section, with and without background 
terms. 
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Fig. 3. Right panel: Pion momentum differential cross section in the LAB frame. Left 
panel: Pion angular differential cross section. 

k/j, and k'^ four momenta are proportional (q 2 = 0) and k n ~ q is assumed 
everywhere except in the nuclear form factor. Thus, the RS differential cross 
section depends on cos 9^ or t only through the nuclear form factor and any 
further cos 6 n or t behaviour induced by the dependence of the amplitudes 
on k n is totally neglected. This is a good approximation at neutrino ener- 
gies above 2 GeV. However, at the energies relevant in the MiniBooNE and 
T2K experiments non parallel configurations become important, and the RS 
model less reliable. We have re-derived RS's expression within our model by 




E^l-cosB^ (GeV) (l-cosS,) (GeV) 



Fig. 4. Laboratory E n (l — cosd n ), at MiniBooNE energies. In the left panel we use 
our full model including full nuclear corrections. In the right panel, we show results 
from the axial contribution of the AP mechanism, neglecting pion distortion and A 
in medium effects. We display the MiniBooNE published histogram (solid), conveniently 
scaled down, from Ref. [1] and MiniBooNE results (dashed histogram) obtained by turning 
off the NUANCE FSI of the outgoing pion (G. Zeller, private communication). 
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considering only the dominant axial part of the AP process (~ C§), neglect- 
ing nuclear medium corrections and replacing k n by q in the pion emission 
vertex. In the right panel of Fig. [4] we see that the new E^fl — cos 9 n ) distri- 
bution is significantly wider than that obtained without implementing this 
replacement and that it reasonably describes the MiniBooNE published dis- 
tribution (solid histogram). The agreement is much better when compared 
with some preliminary MiniBooNE results (dashed histogram) obtained with 
a different treatment of the outgoing pion distortion. This calculation shows 
the uncertainties associated to the t = approximation at low energies. 

Pion distortion induces some additional discrepancies. MiniBooNE im- 
plement this effects through a Monte Carlo cascade model for the tt prop- 
agation in medium. However, coherent cross sections cannot be calculated 
from a Monte Carlo cascade algorithm, because the coherent production is 
a one step process and by using a Monte Carlo algorithm we break the co- 
herence of the process. Nevertheless, one could still reasonably estimate the 
total coherent cross section from the NUANCE FSI cascade if it is used to 
eliminate from the flux not only those pions which get absorbed or suffer 
inelastic processes but also those that undergo QE steps. To our knowledge, 
these latter events are accounted for in the MiniBooNE analysis, despite not 
being coherent. In our calculation the imaginary part of the pion- nucleus 
potential removes from the flux of the outgoing pions those that are absorbed 
or undergo QE interactions. We try to estimate this effect by switching off 
the QE contribution to the pion-nucleus optical potential induced by elastic 
pion-nucleon collisions, and using an optical potential with an imaginary 
part due to absorption and inelastic channels alone. For the MiniBooNE 
flux averaged cross section we find a 20% enhancement (see NC* entry in 
Table [TJ) in good agreement with the effects observed by turning off the NU- 
ANCE FSI. We conclude that the RS model is not as reliable for MiniBooNE 
and T2K experiments as for v energies above 2 GeV. Our model provides an 
E n {l — cos6 n ) distribution much more peaked, and thus it might improve 
the description of the first bin value in Fig. 3b of Ref. [T|. Moreover, the 
drastic change in the £^(1 — cos9 w ) distribution shape might produce some 
mismatch between the absolute normalization of the background, coherent 
and incoherent yields in the MiniBooNE analysis. 

In TableCDwe show our predictions for the MiniBooNE, K2K and T2K [H] 
flux averaged cross sections. Since our model neglects all resonances above 
the A, our predictions become less reliable when the energy increases, so we 
set up a maximum neutrino energy in the flux convolution E msCK , neglecting 
the long tail of the v flux. Up to these energies, one can assume A domi- 
nance and still cover about 90% of the total flux (65% for T2K antineutrino 
flux). We expect corrections (higher cross sections) of around 20-30% to 
our results for MiniBooNE and T2K (larger for K2K). Our prediction lies 
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Reaction 


Exp. 


a[10- 4l W] 


^cxp[10- 4U cm^] 


£ max [MeV] 


CC 




K2K 


4.68 


< 7.7 Q3] 


1.80 


CC 


^+ 12 c 


MiniBooNE 


2.99 




1.45 


CC 


^+ 12 c 


T2K 


2.57 




1.45 


CC 


^+ 16 o 


T2K 


3.03 




1.45 


NC 


^± 12 c 


MiniBooNE 


1.97 


7.7 ±1.6 ±3.6 p] 


1.34 


NC* 


^+ 12 c 


MiniBooNE 


2.38 


7.7 ±1.6 ±3.6 p] 


1.34 


NC 


^+ 12 c 


T2K 


1.82 




1.34 


NC 




T2K 


2.27 




1.35 


CC 


*V± 12 C 


T2K 


2.12 




1.45 


NC 


^+ 12 c 


T2K 


1.50 




1.34 



Table 1. Coherent pion production total cross sections. 



well below the K2K upper bound, while being notably smaller than that 
given in |12j for NC MiniBooNE. However, notice the previous discussion on 
RS model, which is being used in the MiniBoone analysis. The K2K cross 
section and the value quoted in Ref. |l2j seems somehow incompatible with 
the approximate relation acc ~ 2<tnc, expected from A— dominance and 
neglecting finite muon mass effects. 
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